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Abstract. 

A new iterative method for calculating the travel time of a photon as a 
function of the spatial positions of the emitter and the receiver in the field 
of a static, spherically symmetric body is presented. The components of the 
metric are assumed to be expressible in power series in m/r, with m being half 
the Schwarzschild radius of the central body and r a radial coordinate. The 
procedure exclusively works for a light ray the parametric equations of which 
may be expanded in power series of the gravitational constant G. It is shown 
that the expansion of the travel time of a photon along such a ray only involves 
elementary integrals whatever the order of approximation. An expansion of the 
impact parameter in power series of G is also obtained. In order to illustrate the 
method, the perturbation terms are explicitly calculated up to the order G 3 . The 
expressions yielding the third-order terms in the light travel time and the impact 
parameter are new results. 



PACS numbers: 04.20.-q, 04.25.-g, 04.80.Cc, 95.10.Jk 

1. Introduction 

Determining the travel time of a photon as a function of the positions of the emitter 
and the receiver for a given time of emission (or reception) is a crucial problem in 
many tests of general relativity. Indeed, such a function, that is called a time transfer 
function, is relevant for modelling not only the experiments involving the measurement 
of a time delay or the comparison of distant clocks [1,2], but also the bending of light 
and the highly accurate astrometry [3-5]. 

In this paper, we leave out the dynamical aspects which should be eventually 
taken into account in realistic models and we concentrate our attention on the case 
of a static, spherically symmetric metric describing the gravitational held of a central 
body. We suppose that the photon is propagating in a region of space where the 
components of the metric tensor can be assumed to admit analytic expansions in 
powers of m/r, where m is half the Schwarzschild radius of the central body and 
r an isotropic radial coordinate. The metric is thus characterized by an infinity of 
dimensionless constants generalizing the well-known post-Newtonian parameters /3 
and 7. In fact, we restrict our attention to the case where the photon follows a path 
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that we call a quasi-Minkowskian light ray (see [5]), that is a null geodesic described as 
a perturbation in powers of G of a Minkowskian null segment. The corresponding time 
transfer function is then represented by a generalized post-Minkowskian expansion in 
powers of G. For the sake of brevity, a term of order G n will be said to be of order n. 
The first-order term in the expansion of the time transfer function is the well- 
known Shapiro time delay [6], which can be obtained by different reasonings, with 
some of them involving only elementary calculations (see, e.g., [7] or [1]). Nevertheless, 
calculating the higher-order terms is a more difficult problem. As far as we know, two 
distinct approaches are currently available. 

a) The method based on integrating the null geodesic equations. This procedure 
has been applied up to the order G 2 within the framework of the parametrized post- 
post- Newtonian formalism in [8], and then reconsidered in the case of the exterior 
Schwarzschild metric for a class of quasi-Galilean coordinate systems of interest in 
celestial mechanics [9,10]. The problem has been recently revisited in [11] for a three- 
parameter family of static, spherically symmetric spacetimes in the context of the 
Gaia mission. However, the approach developed in [9-11] presents the drawback to 
be indirect, since the results are deduced from a solution which corresponds to a light 
ray emitted at infinity in a given direction. 

b) The methods which are natively adapted to the generic case where both the 
emitter and the receiver of the light rays are located at finite distances from the 
origin of the spatial coordinates. These methods are based either on an iterative 
determination of the Synge or bi-scalar world function (see [12] for the Schwarzschild 
metric and [3] for a more general and simpler approach), or on an iterative integration 
of the Hamilton- Jacobi equation satisfied by the time transfer function (see ref. [13]). 
The case n — 2 has been successfully solved using the two variants. 

The calculations required by the procedures mentioned in a) and b) get quickly 
unwieldy as n is increasing. For this reason, we develop here a new method enabling 
a systematic determination of the time transfer function up to any given order n. 
This approach relies on the null geodesic equations written in spherical variables. An 
integral expression of the photon travel time involving the impact parameter of the 
light ray can be obtained when the photon is not passing through a pericentre. This 
expression transforms into an intcgro-diffcrential equation when the property of the 
impact parameter to be a derivative of the time transfer function is taken into account 
(see [5]). We show that this equation can be solved by an iterative procedure and that 
each perturbation term involved in the time transfer function is expressible as a sum of 
elementary integrals which are easy to calculate with any symbolic computer program. 
A theorem of analyticity of the perturbation terms proved in the present paper allows 
to extend the results to the more general cases where the light ray is passing through 
a pericentre. Since the property of analyticity we use is inferred from the above- 
mentioned Hamilton- Jacobi equation, it may be pointed out that our procedure is an 
hybrid of the above-mentioned methods. 

In order to illustrate the convenience of the new method, we carry out the explicit 
calculation of the time transfer function up to the order G 3 . The formula obtained here 
for the third-order term is new. It must be emphasized that our formula substantially 
improves the results previously obtained in [14] and [15]. Indeed, the expressions 
of the light travel time obtained in these works involve the radial coordinate of the 
pericentre of the ray without explicitly calculating this quantity as a function of the 
spatial positions of the emitter and the receiver. 

It may be argued that the most accurate projects for testing general relativity, like 
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LATOR [16] or ASTROD [17], are generally considered as only requiring a knowledge 
of the propagation of light up to the second order (see, e.g., [18] and references 
therein). However, the possibility for an effect of order n + 1 to be greater than some 
contributions of order n cannot be excluded, as it is pointed out in [19]. The necessity 
to be cautious has been recently shown for the Gaia mission. This mission, indeed, 
is currently tackled within the linearized, post-Minkowskian regime [20, 21] or the 
usual post-Newtonian approximation [22,23]. Nevertheless, an apparent discrepancy 
between the standard approach and the numerical integration of the null geodesic 
equations has required an in-depth discussion of an 'enhanced term' of order G 2 
(see [11] and [5]). For this reason, it seems to us that exploring new systematic 
procedures enabling the calculation of the time transfer function at any order is fully 
justified. 

The paper is organized as follows. Section 2 lists the notations and conventions 
we use. In section 3, the general assumptions on the metric are stated. Section 4 
is devoted to the expansion in a series in powers of G of the time transfer function 
associated with a quasi-Minkowskian light ray. Section 5 yields a recurrence relation 
satisfied by the perturbations terms involved in the expansion of the time transfer 
function. A fundamental property of analyticity is established for these terms in 
section 6. The impact parameter of any non-radial quasi-Minkowskian light ray is 
shown to be expressible as a series in powers of G in section 7. This feature is the basis 
of the new iterative procedure proposed in this paper for determining the time transfer 
function at any order n. This procedure is implemented in section 8. A simplification 
is carried out in section 9. The expression of the time transfer function is given up 
to the third-order in G in section 10. The impact parameter of the corresponding 
light ray is obtained up to the same level of approximation in section 11. The impact 
parameter of a ray emitted at infinity is obtained in section 12. Concluding remarks 
are given in section 13. An appendix yields some hints for the hand calculation of the 
third-order term in the time transfer function. 

2. General assumptions, notations and conventions 

Our general assumptions, notations and conventions are the following. 

• Spacetime is assumed to be a static, spherically symmetric manifold (Vi,g). 
We suppose that there exists a region T>h in which the metric g is regular, 
asymptotically flat and may be interpreted as the gravitational field of a central 
body having a mass M. We put m = GM/c 2 . 

• We assume that T>h may be entirely covered by a single quasi-Cartesian system 
of coordinates (x°,x l ) adapted to the symmetries of the metric. We put x° = ct, 
with t being a time coordinate, and x = (x l ). 

• Greek indices run from to 3, and latin indices run from 1 to 3. 

• The signature adopted for the metric is H . 

• Any bold italic letter refers to an (ordered) triple: (a 1 , a 2 , a 3 ) = (a 1 ) = a. All the 
triples are regarded as 3- vectors of the ordinary Euclidean space. 

• Given two triples a and b, a.b denotes the Euclidean scalar product a l b % , with 
Einstein's convention on repeated indices being used. 

• \a\ denotes the formal Euclidean norm of the triple a: \a\ — (a. a) 1 / 2 . If \a\ = 1, 
a is conventionally called a unit (Euclidean) 3- vector. 
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• a x b is the triple obtained by the usual rule giving the exterior product of two 
vectors of the Euclidean space. 

• Given a bi-scalar function F(x,y), X? x F(x,y) and 'V y F(x,y) denote the 
gradients of F with respect to x and y, respectively. 

3. Generalized post-Minkowskian expansion of the metric 

For convenience, the coordinates (x , x) are chosen so that the metric takes an isotropic 
form in the domain of regularity T>h : 

ds 2 = A(r)(dx ) 2 - B- 1 (r)S iJ dx i dx : ', (1) 

where r = \x\. Using the corresponding spherical coordinates (r, $, ip), one has 

S ij dx i dx j = dr 2 + r 2 dd 2 + r 2 sin 2 Mip 2 . 

The light rays of the metric (1) are also the light rays of any metric ds 2 conformal 
to (1). This feature enables us to carry out our calculations for a metric containing 
only one potential. We choose ds 2 = A~ l (r)ds 2 , that is 

dS 2 = (dx°) 2 -U(r)8 l2 dx l dx\ (2) 

where IA is defined by 

u{r) = WWy (3) 

The metric (1) is considered as a generalization of the exterior Schwarzschild 
metric, which may be written in the form 



1 



it) 



2 



ds 2 = r 2rJ {dx 0)2 _ f x m\ s dx i dx j (4) 

/ m\ 2 V 2rJ 

\ + 2r) 
in the region outside the event horizon located at r = to/2. So we henceforth admit 
that there exists a value rh > of the radial coordinate such that the domain of 
regularity T>u is the region outside the sphere of radius rh- If the metric admits at 
least one event horizon, we must take for rh the value of r on the outer horizon. By 
analogy with general relativity we consider that rh ~ to and we assume that whatever 
r > rh, A(r) and B~ l (r) are positive functions represented by analytical expansions 
as follow: 

., . _, 2to „„to 2 3 . m 3 to 4 v^ (— l)"n to" ,_, . 

n— 5 
o 2 1 3 1 4 °° n 

„ i. . m 6 m 1 m 1 m v—v, x m , ,. 

B_1 W = ! + 2 ^7 + 2 £ — + 2 73 ~ + 17^— + E(> - !) — - ( 5& ) 

n— 5 

where the coefficients /3, /3 3 , . . . , /3 n , 7, e, 73, . . . , 7„, . . . are generalized post-Newtonian 
parameters chosen so that 

/3 = 7 = e = 1, ^„=7„ = 1 for n>3 (6) 

in general relativity. It results from (3), (5a) and (56) that the potential U(r) occurring 
in (2) may be written as 

00 n 

U(r) = 1 + 2(1 + 7)^ + ^2^^7, (7) 

n— 2 
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for r > r/j, with the coefficients K n being constants which can be expressed in terms 
of the generalized post- Newtonian parameters involved in the expansions of A(r) and 
B[r). Taking into account a notation already introduced in [5], namely 

K = 2(l + 7)-0+!c, (8) 

K2 and K3 are given by 

k 2 = k, k 3 = 2k- 2/3(1 + 7 ) + |^3 + I 73 . (9) 

In general relativity, we have 

K 2 = K=f, K 3 = §. (10) 

4. Time transfer function for a quasi-Minkowskian light ray 

Let x A and x B be two points located in T>h- Throughout this paper, we generically 
consider a photon emitted at x A and received at x B . The time of emission and the 
time of reception of this photon are denoted by t A and t B , respectively. It is assumed 
that the light ray followed by the photon is a null geodesic path T which is entirely 
lying in the domain of regularity V^. As it has been emphasized in introduction, it 
would be of primordial interest for modelling a lot of relativistic effects to determine 
the time transfer function associated with T, that is, the expression giving the travel 
time t B — t A as a function of x A and x B : 

t B -t A =Tb(x A ,x B ). (11) 

In practice, however, the problem is extremely complicated since there exists in 
general an infinite set of light rays emitted at x A at a given time t A and passing 
through x B (see, e.g., [24,25] for the Schwarzschild metric and [26] for a larger class 
of spacetimes). So, in this paper, we restrict our attention to the special class of null 
geodesic paths we have called the quasi-Minkowskian light rays in [5] . This means that 
in what follows, the path covered by the photon is assumed to be entirely confined in 
T>h and to admit parametric equations of the form 



x° = ct A +C\x B - x A \ +^X'( n) {x Al x B ,C) 1 (12) 

71=1 

OO 

x =z(0 + ^2x (n) (x A ,x B ,£), (13) 

71=1 

where £ is the affine parameter varying on the range < £ < 1, z(£) is defined by 

z(£) =x A +£(x B -x A ) (14) 

and the functions X? -, and -XVn) are terms of order n obeying the boundary conditions 

X ( °„ ) (x A ,x B ,0)=0, (15) 

X( n )(x A ,x B ,0) = X {n) (x A ,x B ,l) = 0. (16) 

According to a notation already introduced in [5], such a null geodesic path will be 
denoted by T s (x A , x B )\. For the sake of brevity, the time transfer function associated 
with T s {x Al x B ) will be henceforth denoted by T{x Al x B ) or simply by T. Setting 

\ In a static spacetime, the mention of the initial time t A may be omitted. 
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£ = 1 in (12), it may be seen that this function admits an expansion in power series 
of G as follows: 



T(x A ,x B ) = lXB Xa] +T / r {n \x A ,x B ), (17) 



where T^ n > stands for the term of order n. 

Expansion (17) is easy to determine when x A and x B are linked by a radial null 
geodesic entirely lying in V^. In this case, indeed, it is immediately deduced from (2) 
that the expression of T is given by the exact formula 

1 



T(r A ,r B ) = sgn(r B - r A )- / y/U{r)dr, (18) 



C -'r A 



where r A = \x A \ and r B = \x B \. Substituting for U(r) from (7) into (18) shows that 
T may be expanded as follows: 



oc 



T(r,,r B ) = ^^ + VT(")(r A ,r B ), (19) 



c £i 



where the first three perturbation terms are given by 

T^\r A ,r B )= { ±±J^ 



ln^ 

r A 



(20) 



T^(r A ,r B )=[ K -l(l + ~ f f]^-^-^, (21) 

r( 3 )(r„r B ) = I[ K3 -(l+7)« + i(l+7) 3 ]^f- + -)^^^. (22) 

r A r B \r A r B J c 

Determining the right-hand side of (17) is much more complicated when 
Y a (x Al x B ) is not a radial geodesic. As it has been recalled in introduction, the 
perturbations terms 7"'™' might be obtained by an iterative integration of the null 
geodesic equations. Indeed, taking into account that ds 2 = along a null geodesic, it 
results from (2) and (13) that the time transfer function is given by 



1 /-i , 

T(x A ,x B ) = - / VWW) 

cJo 



EdX( n ^(x A , x B , £) 
dl 



d£, (23) 



where the integral is taken along T s (x A ,x B ). Taking into account the boundary 
conditions (16), it may be inferred from (23) that each function T*-"-* is theoretically 
calculable if the perturbations terms -X^i), . . . , X/ n _ 1 - ) involved in (13) are determined 
by solving the null geodesic equations. This procedure is cumbersome, however. So 
we develop a different method, based on a property of analyticity of the functions J~y n > 
which can be deduced from a recurrence relation (see sections 5 and 6). 

5. Recurrence relation satisfied by the functions T*-"^ 

It follows from a general result established in [13] that T(x A ,x) satisfies a Hamilton- 
Jacobi equation which is equivalent to the eikonal equation 

c 2 \V x T(x A ,x)\ 2 =U(r) (24) 

when the metric is given by (2). Replacing T(x A , x) in (24) by its expansion in powers 
of G and U(r) by (7), and then applying the formulae already found in [13], we get a 
proposition as follows. 
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Proposition 1 The perturbation terms T^ n > involved in expansion (17) may be 
written in the form 

T {n \x A ,x B ) = -\x B - x A \F^(x A ,x B ), (25) 

c 

where the functions F^ n > are determined by the recurrence relation 

F^{x Al x B ) = {l + 1 )m / i-^rr, (26) 

Jo \ z (t,)\ 



fW(x ^^ K " ffl "(kiF 



„2 n ~ x r 1 



r 2 f r i 

-tE / V x T {p) (x A ,x).V x T {n - p) {x A ,x) dH (27) 

/or n > 2, tuiift z(£) feeing defined by (14)%- 

The recurrence relation explicited in proposition 1 shows that T s (a; A ,a; B ) is 
unique provided that expansion (17) is an admissible representation of the time 
transfer function. However, determining the most general conditions under which 
our construction is valid remains an open problem. According to [13], (26) and (27) 
are inferred from an integro-differcntial equation involving the analytic expansion of 
the metric along the straight segment joining x A and x B . Consequently, we shall 
henceforth assume that the expression \x B —x A \{l + ^2™ 1 F^ p \x A ,x B )]/c constitutes 
a reliable approximation of the time transfer function as long as the straight segment 
joining x A and x B does not intersect the hypersurfacc r = r>,, a condition expressed 
by the inequality 

\z(0\>r h for 0<<£<1. (28) 

This condition is largely satisfied for a star observed in the Solar System, as well as 
in any foreseeable test of general relativity in the vicinity of the Sun. 

It has been shown in [13] that a recurrent relation equivalent to the one stated by 
proposition 1 enables to carry out the calculation of T*- 1 -* and 1~( 2 \ As matters stand, 
however, we do not know whether (27) allows explicit computations for n > 3. Going 
deeper into this question is beyond the scope of this paper. Proposition 1 is used here 
only for proving a property of analyticity which is indispensable for justifying our new 
procedure. 

6. Analyticity of the functions T^ 

Let us begin with proving the following lemma. 

Lemma 1 The functions F^ n > recursively determined by (26) et (27) are analytic in 
x A and x B , except when x A and x B are such that n B — —n A , with n A and n B being 
defined as 

n A = — , n B = — . (29) 

r A r B 

§ Let us emphasize that the integrals involved in (26) and (27) are taken along the straight segment 
admitting the parametric equation x = z(£), < £ < 1. 
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Proof of lemma 1. The proposition is obviously true for n = 1, since the 
integrand l/|z(£)| in (26) is analytic in x A and x B for any £ such that < £ < 1, 
provided that n B ^ — n A . Suppose now the validity of lemma 1 for F^\ . . . ,F^ n \ 
Assuming p to be such that 1 < p < n, and then substituting z(£) for x into 
^xT~( p \x A , x), it is immediately inferred from (25) that 



V x T {p Hx A ,x) 



Jx=z(Q 

+£,\x B - x 
where N AB is defined by 



= N ab fM(x a ,z(£)) 

V x F^\x A ,x 



«(« 



N AB = 



*£ B **s J- 

\X R X l 



(30) 
(31) 



Using (30) leads to 

,2 



F<")(a:x ) aOV 1B F< n+1 -i , >(a: A ,a:) 



v x r (p) (^^)-v a! r ( " +1 - p) (^^) 

+C(a; B -a; 

+ F(" +1 -p)(x j4 ,x)V xj fW(x a ,x) 



= F^{x A ,z(0)F (n+1 -vHx A ,z(0) 



x=z(0 

V x F^(x A ,x).V x F^ n+1 -P\x A ,x) 



*(«) 



(32) 



It follows from our assumption that the right-hand side of (32) is a sum of functions 
which are analytic in x A and x B for any £ such that < £ < 1, except if n B = — n A - 
Each integral 



V x T (p) (as,. , x). V x T {n+1 - p) (x A ,x] 



x=z (i) 



dt 



is therefore analytic if n B ^= —n A . The same property is obviously possessed by the 
integral L d£/|z(£)|" +1 . Lemma 1 is thus proved by recurrence. 

Since \x B — x A \ is analytic except if x B =/= x Al we can state the proposition below. 

Proposition 2 The functions 7"'™' involved in expansion (17) are analytic in x A and 
x B when both the following conditions are met: 

a J x b 7^ X A ; 

b) n B ^ -n A . 

The importance of this property will clearly appear in section 8. It is worth 
of noting that the second condition in proposition 2 is automatically fulfilled when 
inequality (28) is satisfied. This fact explains why the condition b) is never explicitly 
involved in the assumptions of the propositions enunciated below. 

7. Relation between T(x A ,x B ) and the impact parameter of T s (x A , x B ) 

7.1. Null geodesic equations 

Let r be an arbitrary non-radial null geodesic path of the metric ds 2 . We suppose 
that r is confined in the region T>h and described by parametric equations x a = x a ((), 
where £ is an arbitrarily chosen affine parameter. We choose the spherical coordinates 
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(r, i?, ip) so that $ = 7r/2 for any point of this path. Denoting by l a the covariant 
components of the vector tangent to r,j(a; A , x B ), an equation as follows 

lodx + l r dr + l v dip = (33) 

is satisfied along T since l a is a null vector. Owing to the symmetries of the metric, 
we have 

h = E, (34) 

l v = -J, (35) 

with E and J being constants of the motion. For convenience, the affine parameter 
is chosen in such a way that E > 0. Furthermore, it is always possible to suppose 
J > without lack of generality when calculating the time transfer function in a 
static, spherically symmetric spacetime. Then the quantity defined as 

*>=| (36) 

is the impact parameter of the light ray (see, e.g., [27] and [5])||. It may be noted that 
6 = would correspond to a radial null geodesic. 

Since ds 2 = along T, it follows from (34), (35) and (36) that 

E 



l r = -e-Jr 2 U(r)-b 2 , (37) 

r 

where e = 1 when r is an increasing function of time and e = — 1 when r is a 

decreasing function of time^f. Substituting for l r from (37) into (33), and then dividing 

throughout by E, we get a relation enabling to determine the time transfer function 

by an integration along r s (a;_4,a; B ), namely 



dx° = bdip+- Jr 2 U(r) - b 2 dr. (38) 

r 

Our procedure for calculating explicitly each function 7"'™' rests on the property 
shown in the next subsection that the impact parameter b can be determined as a 
function of x A and x B by taking the partial derivative of T with respect to the cosine 
of the angle formed by x A and x B . 

7.2. Expansion of the impact parameter as a series in powers of G 

Let [ip A , ip B ] be the range of the angular function ip(t) along a quasi-Minkowskian light 
ray T s (x A , x B ). For the sake of brevity, we shall frequently use a notation as follows 

/i = n A .n B = cos(<,5 B - ip A ). (39) 

Using this notation, the time transfer function may be considered as a function of 
r A ,r B and ^: 

T(x A , x B ) = T(r A , r B ,n). 

It is then possible to enunciate the following proposition. 

|| b is an intrinsic quantity attached to T since the constants of the motion E and J are themselves 
coordinate-independent quantities. 

% The sign of e in equation (37) is changed if and only if the photon passes through a pericentre or 
an apocentre. The passage through an apoccntre corresponds to an extreme relativistic case. 
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Proposition 3 Let x A and x B be two points in T>h such that both the conditions 
n A 7^ n B an d (28) are fulfilled. The impact parameter b of a quasi- Minkowskian light 
ray joining x A and x B may be expanded in powers of G as follows: 



i+E 



(40) 



where r c is the usual Euclidean distance between the origin O of the spatial coordinates 
and the straight line passing through x A and x B , namely 

r c = -J , \n A x n B \, (41) 

\x B x A \ 

and the quantities q n are functions of x A and x B given by 



q n = -c ( r -l) n VEZ aTMfrx.r.,/,) 

Ira/ r c d/j, 

Proof of proposition 3. Noting that 
\n A x n B \ = \J\ - n 2 , 

it is immediately inferred from equation (13) in [5] that the impact parameter of 
T(x A , x B ) may be rewritten in the form 

■dT(r Al r B ,fi) 



b = -c\J\- ix z -z • (43) 

Substituting for T from (17) into (43) directly leads to the expansion given by (40). 
The zeroth-order term is easily derived from the elementary formula 



\x B -x A \ = \jr\ - 2r A r B [i + r\. (44) 

Indeed, using this expression and taking (41) into account yield 

a\x B — x A \ _ r c . . 

We shall see in the next section that the expression of the time transfer function 
corresponding to a quasi-Minkowskian light ray can be straightforwardly deduced from 
proposition 3. 

8. Implementation of the method 

If r a (x A ,x B ) passes through a pericentre x Pl the integration of (38) requires the 
determination of the value of the radial variable at x P as a function of x A and x B . The 
calculation of the time transfer function is very complicated for such a configuration. 
Fortunately, owing to the analytic extension theorem, it follows from proposition 2 
that it is sufficient to determine the expression of each term T( n > as a function of x A 
and x B in an arbitrarily chosen open subset of the domain of analyticity. For this 
reason, the calculation of the T^ are henceforth carried out under the assumption 
that x A and x B fulfil the following conditions: 

a) The radial variable r along a quasi-Minkowskian null geodesic joining x A and 
x B is an increasing function of t: 

dr 

— >0, t A <t<t B . (46) 
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b) An inequality as follows 

N AB .n A > (47) 

is satisfied, with N AB being denned by (31). 

These conditions considerably simplify the calculations. Indeed, (46) eliminates 
the occurrence of any pericentre (or apocentre) between the emission and the reception 
of light and (47) implies that the projection of the origin O on the straight line passing 
through x A and x B lies outside the straight segment linking x A and x B . One has 
therefore 

r c < r A < r < r B (48) 

for any point of F s (x A , x B ). These inequalities ensure that condition (28) is met, since 
r h < r A for any point x A located in 2\. 

Under these assumptions, integrating (38) along T s (x A ,x B ) is straightforward 
since the range of the angular function <p(r) between the emission and the reception 
of the photon is given by 

<Pb — <Pa = arccos \x. (49) 

Noting that in this case s = 1, it may be seen that the time transfer function is then 
related to the impact parameter b by an equation as follows 



/ (x Al x B ) — — 
c 



r r B i ■ 

b arccos /J.+ — \Jr 2 U(r) — b 2 dr 

J r a r 



(50) 



Since b is a function of x A and x B determined by (43), (50) has to be regarded 
as an integro-diffcrential equation satisfied by T. In order to solve this integro- 
diffcrential equation by an iterative procedure, let us substitute (7) for U and (40) 
for b. Expanding yfr 2 U(r) — b 2 jr in a power series in m/r c , rearranging the terms 
and introducing the notation 

s = yr 2 — r 2 , (51) 

we get an expression as follows for T 



I ( X Al X BJ — 

c 



r c arccos /i + / - dr 




1 

r c q n arccos /i 

' r A 



U n - 



rhn 



dr}, (52) 



where each U n is a function of r which may be written in the form 

K = £±&1, (53) 

s 

3 «-4 h-n+1 

U„=J2 U ^(Qu ■ • •,? n -i)r* l - fc - 2 -2 S=r (54) 

fe=0 

for n > 2, with the quantities Ukn{qi, ■ ■ ■ ,<Zn-i) being polynomials in qi,...,q n —i. 
Noting that 

fTB 

r c arccos /i + / - dr = \x B — x A \ (55) 

and 

arccos [x — r c I — = (56) 



+ 


/ - dr = \x 
>r A r 


r c 


lr A rs 
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when conditions (46) and (47) are met + , (17) is immediately recovered from (52), with 

each perturbation term being given by 

1 / \ n f rB 
T^(x A ,x B ) = -(-) / U n dr. (57) 

c \ r cJ J rA 

Equation (57) constitutes the main ingredient of the procedure developed in this 
paper. The fact that the coefficient q n is not involved in U n and the property for 
each coefficient q^ to be proportional to a derivative of the function T^> imply that 
T^ n > can be determined when the sequence of functions T^ l \ . . . , X'™ -1 -' is known. 
Moreover, it follows from (54) that all the integrations involved in the right-hand side 
of (57) are elementary and can be carried out with any symbolic computer program. 
Consequently, our procedure enables us to perform the explicit calculation of 7"'™' as 
a function of x A and x B whatever the order n. 

The problem is treated here in a detailed manner up the third order. So, (53) 
must be supplemented by the expressions of U 2 and U 3l namely 

(58) 



u 2 = 


4 
rs 3 


u 3 = 


K 3 r 7 c 
r 2 s 5 



(1 + T)< 


h r 3 t [2k- 


(1+7) 2 - 


q\\r\r 


s 3 


[2k 3 - (1 + 


2s 3 
7)(« + 9i " 


- Q2M 


rs 5 
[2 K 


-3(l + 7 ) 2 


s 5 
- <ll + 2q 2 ] 


qir^r 



2s 5 
[2k 3 - (1 + 7 )(2/s - q\ - 2q 2 ) + (1 + ~/) 3 }r 3 c r 2 qi q 2 r 2 r 3 



(59) 



2s 5 s 5 

The expression of T' 1 ' is immediately inferred from (53) and (57). Noting that 



\Z r A - r c = r A N AB .n A , (60a) 

\Jr\ - r 2 = r B N AB .n B (606) 

when conditions (46) and (47) are met, we get 

n-iDf \ (1+7)™, r B +N AB .x B 

T () (x A7 x B ) = In — — — . (61) 

c r A + N AB .x A 

As it could be expected, we recover the well-known Shapiro time delay expressed in 

isotropic coordinates (see, e.g., [7]). 

For determining q\, it is preferable to rewrite (61) in the more elegant form (see, 

e.g., [1]) 

T W( Xa ,x b ) = ii±^ln ( r *+ r ° + \*s-*A . (62) 

c \r A +r B - \x B -x A \) 

Substituting for T^ 1 ' from (62) into (42) written for n — 1, and then using (45), it is 
easily seen that 

1 + n A .n B \r A r B J 
Taking into account this determination of gi, it would be possible to carry out the 
calculation of T^ via (58). Then, q 2 could be derived from (42) taken for n = 2. 
Consequently, T*- 3 -* could be deduced from (59). However, we shall see in the next 
section that the method can be simplified by making use of the differential equation 
governing the variation of the angular coordinate along the light ray. In particular, it 
turns out that determining q 2 is not indispensable for calculating T^ 3 > . 

+ Note that (55) is just (50) written in the case where the gravitational field vanishes, i.e., m = 0. 
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9. Simplification of the procedure by using a constraint equation 



Equations (35) and (37) are equivalent to the geodesic equations 
dip J 

~d( = r 2 U{r) ' 
dr E 



(64) 



dC -rUir)^^- (65) 

Eliminating the affinc parameter £ between (64) and (65) leads to 

d ^=e b - - 1 (66) 

dr r ^Jr 2 U{r) - b 2 

Since s = 1 when conditions (46) and (47) are met, integrating (66) and taking 
into account (49) yield thereby 

rr B b 

arccosu= / dr (67) 

J rA r^Jr 2 U{r) ~ b 2 

Equation (67) implicitly determines b as a function of x A and x B . So it may be 
expected that this equation implies some conditions on the coefficients q n which may 
be used to simplify the calculations. 

Replacing U by (7) and b by (40) into (67), it may be seen that 



dr 
arccos /i 



r A TS ^S^' 



T E 7 fWndr, (68) 

where the W m 's are functions of r which may be written in the form 

(69) 

9 
TV 

>-x ( 7 °) 



W-i. = 


3 2 

= -(1+7)% + 91^, 




3(n-l) 


w n -- 


= Yl w kn(qi,---,q n -i)rl n 




fe=0 



„fc— n+1 

„2n+l 



for n > 2, with the terms Wfc n (<Zi, • • • , Qn-i) being polynomials in qi, . . . , g n _i. Taking 
into account (56), it is immediately seen that (68) reduces to 



m x ,TB 



E r / ^ dr = ( 71 ) 



n=l \ rc/ ""-a 

for n > 1. Since (71) holds whatever m, it is clear that (68) is equivalent to the infinite 
set of equations 

W n dr = 0, n = l,2,... (72) 

The set of constraint equations (72) may be systematically used for simplifying 
our problem. Let us consider the functions U* defined as 

Ul = U u (73) 

ra-l 
p=l 
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for n > 2, where the fc p „'s are arbitrary quantities which do not depend on r. Taking 
into account (72), it is immediately seen that 

U n dr= / U*dr. (75) 

ir A Jr A 

Hence T^ n > may be rewritten in the form 

T^(x A ,x B ) = -(y) I U* n dr. (76) 

It is easily seen that a judicious choice of the quantities k pn enables to shorten 
the expressions involved in (76) when n > 2. Until n = 3, only the expression of W\ 
is needed. Indeed, for n = 2, setting fci 2 = 5(71 removes the term in q 2 and leads to 

r » _ « r c , (l + 7)girg , [2^-(l+ 7 ) 2 ]r 2 r 

^ 2 " 7^ + 2i3 + 2^ • V 1 ' 

For n = 3, choosing fci3 = (72 and fc 2 3 = remove the terms involving q 2 . Then U£ 
reduces to 

K 3 r 7 c nqirl [2k 3 - (1 + j)(k + ql)]r 5 c 



US = 



[2«-3(l+ 7 ) 2 -g?]gir*r 

1 -I- ^/) 3 ]r 3 r 2 

(78) 



2s 5 
[2k 3 - (1 + 7) (2k - g 2 ) + (1 + -ff}r^r 2 



2s 

It is thus proved that owing to the constraint equation (67), only the 
determination of q\ is required for calculating the functions T^ 2 -* and T^ ■ 

Remark. It may be pointed out that the coefficients q n could be directly inferred 
from the constraint equation without differentiating the functions T^ n > with respect 
to [i. Indeed, it follows from (69), (70) and (72) that 



l+7^y/r| -r 2 -r B y/r 2 - r 2 
1 v/r 2 - r 2 ^Jr\ - r 2 



r c y/r 2 -r 2 - ^/r 2 - r 2 

3 ( n - 1 ) rr B k-n+1 

X J2 W kn(Qu---,Qn-i)rl n - k - 1 / ~^rdr (80) 

fc=0 Jr A s 

for n > 2. Equation (80) shows that q n can be determined once Qi, • • • ,q n —i arG 
known. 

It is easily checked that (79) is equivalent to (63). Indeed, noting that 

V r l -r 2 c - \jr 2 A -r 2 = \x B - x A \, (81) 

and then taking into account (41), (60a) and (60&), it may be seen that (79) transforms 

into 

N AB .n B -N AB .n A 

<7i= 1+7) i i • 82) 

\n A x n B \ 

Substituting r A n A for x A and r B n B for x B into the numerator of the right-handsidc 
of (31) yields 

\r „ nr „ (r A +r B )(l-n A .n B ) . . 

N AB .n B - N AB .n A = : : . (83) 

\X R X A 
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Finally, substituting for N AB .n B — N AB .n A from (83) into (82), and then noting that 
(41) is equivalent to 

1 r, 1 



\x B -x A \ r A r B \n A x n B \ 
it is immediately seen that (63) is recovered. 

10. Time transfer function up to the third order 

We are now in a position to determine the perturbation terms involved in the expansion 
of the time transfer function up to the order G 3 . The term T^ 1 ' has been already 
treated in section 8. For n — 2 and n — 3, it follows from (77) and (78) that T^™-* may 
be written in the form 

1 /m\ n ' 7 J~} r rB r k - n + 1 

T{n) ^^) = \{^) E U *Mrl n - k - 2 l A ~^rdr (84) 

where <r(2) = 2 and cr(3) = 4, with the coefficients [/£ being polynomials in q\. 
The integrals occurring into the right-hand side of (84) arc elementary and can be 
expressed in terms of r A , r B , r c , \Jr\ — r 2 and \/r B — rj? . For the explicit calculations, 
it is convenient to write (60 a) and (606) in the form 

r~a n r A [r B n — r A ) . . 

^^= \ Xb - Xa \ ■ ( 85fl ) 

r B (r B -r A [i) 



V^^! = ,;;_, J - (856) 

Using (41), (63), (85a) and (856), it may be seen that T*- 2 -* and T^ can be expressed 
in terms of r A r Bl l/r A + l/r B , \x B — x A \ and /i. It has been already emphasized 
that the explicit calculations can be performed with any symbolic computer program. 
Of course, a simple hand calculation is also possible. For n = 2, the result is 
straightforwardly obtained. For n = 3, however, the calculations are somewhat lengthy 
and tedious. For this reason, some hints concerning this case are delivered in the 
appendix. 

We have seen in section 8 that the expressions thus obtained can be considered 
as valid even when conditions (46) and (47) are not fulfilled. So we can formulate the 
following proposition. 

Proposition 4 Let x A and x B be two points in T>h such that both the conditions 
n A 7^ n B and (28) are met. Then T^ 2 ' and T 3 ' are given by 

axccosn A .n B (I+7) 2 



/ (x Al x B ) 



r A r B c 



K- 



<t-(3) / \ m ( l l \ \x B x 

I (x A ,x B ) 



n A x n B \ 1 + n A .n 

K3 - (l + 7)« 



(86) 
arccos n A .n B 



r A r B \r A r B ) c(l + n A .n B ) 

(1+7) 3 
1 + n A .n 

where the coefficients k and K3 are determined by (8) and (9), respectively. 



n A x n B \ 

(87) 
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The formula obtained for T*- 3 -* is new*. Concerning T*- 2 -*, it appears that (86) 
coincides with the expression previously obtained by completely different methods 
in [3] and [13] (j. This concordance confirms the reliability of the procedure presented 
here. 

In the case of a radial null geodesic, it is immediately inferred from 

arccosn A .n B 

hm — | 1— = 1 (88) 

n B ^n A \n A x n B \ 

that the system (20)- (22) is recovered from (62), (86) and (87). This agreement 
constitutes another confirmation of the validity of the new procedure. 

In general relativity, the expressions of T^\ T^ and T*- 3 -* are obtained by setting 
7 = 1, K = ^ and n 3 = |. 

11. Impact parameter up to the third order 

The coefficient q\ has been previously inferred from (62) in section 8. Substituting 
for T^ 2 -* from (86) into (42), and then taking into account that r c \x B — x A \/r A r B = 



(89) 
(90) 



y/1 -^ 2 , we get 




n,-„ l L arccos ^ n , (!+7) 2 n „. ,2 n A 


h l^-^PiVl-ZJ 2 ! + M U ' J j' 


where Q AB is defined as 


Qab = \x B - x A \ 2 n - r A r B (l - n 2 ). 


Then, substituting for T^ 


from (87) into (42) and using (90), we obtain 

, 3 [\x B -x A \ 2 - Q AB ] 


_ r A +r B y/1 - fi J^ 

\X B -X A \ 3 y/l + H \ 


-(l + 7)« 


, ,2 , \ x b — X A \ [1 — (J,) — Qab 


l ~ V±-v 2 J 


+ (l+7) 3 ' 


x B - a; A | 2 (2 - n) - Q AB \ 



( 91 ) 
l + /i J x ' 

The dimensionless coefficients q n can be expressed in terms of the sine (or cosine) 
of the angles formed by n A , n B and N AB . Noting that 

r c (— + —)= \N AB x n A \ + \N AB x n B \ (92) 

\r A r B J 

and that (90) may be written in the form 

Qab = \x B - x A \ 2 (N AB .n A )(N AB .n B ), (93) 

a proposition as follows is straightforwardly inferred from (63), (89) and (91). 

* We have given the expression of 7~' 3 ' without demonstration in [28] . 

J See also [11], where the expression of T^ 2 ' is obtained in harmonic coordinates by an integration 

of the null geodesic equations. 
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Proposition 5 Under the assumption of proposition 4, the coefficients q\, q 2 and q 3 
involved in the expansion of the impact parameter b of a quasi- Minkowskian light ray 
joining x A and x B are given by 

i \ m , ,\N AB x n A \ + \N AB x n B \ 

qi{x A ,x B ) = (1+7) J rv ' ( 94 ) 

1 + n A .n B 



q 2 (x A ,x i 



(N AB .n A )(N AB .n B ) 
1 : : arccos n A .n B 



2 1 - (N AB .n A )(N AB .n B ) 



-(1 + 7) 2 v T — "— ' , (95) 

1 + n A .n B 

g 3 (a; A ,a; B ) = — ^ k 3 1 - (N AB .n A )(N AB .n B )\ 

l + n A .n B 



-(l + 7)« 



1 - n A .n B - (N AB .n A )(N AB .n B ) 



arccos n A .n B 



3 2-n A .n B - (N AB .n A )(N AB .n t 



l + n A .n B J 

Equations (94) and (95) arc identical to the expressions obtained in [5]. We 
remark that the first-order expression of b yielded by (40) and (94) coincides with the 
Euclidean norm of the vector d' given by equation (62) in [11]. On the other hand, 
the formula (96) is new and completes the implementation of our method up to the 
third order. 

12. Case of a light ray emitted at infinity 

A quasi-Minkowskian light ray coming from infinity in a initial direction defined by 
a given unit vector N e and observed at a given point x B is a relevant limiting case 
for modelling a lot of astrometric measurements. According to a notation introduced 
in [5], such a ray is denoted by T s (N e , x B ). The corresponding null geodesic is assumed 
to be a perturbation in powers of G of the straight segment defined by the parametric 
equations 

x° {0 )(X) ^ ct B + Xr c , X( )(A) = Ar c 7V e + x B , -oo < A < 0, (97) 

where 

r c = r B \N e xn B \. (98) 

In order to ensure that condition (28) is satisfied for any point x A of T s (N e , x B ), 
N e and x B must be supposed to satisfy the condition 

|Ar c iV e + x B \ > r h when - oo < A < 0. (99) 

This condition means that the straight segment coming from infinity in the direction 
N e and ending at x B is entirely lying in T>h- 

Then the following proposition can be stated. 

Proposition 6 Let N e be a unit vector and x B a point in T>h fulfilling condition 
(99). The impact parameter of a quasi-Minkowskian light ray emitted at infinity in 



New method for determining the light travel time... 18 

the direction N e and arriving at x B is given by expansion (40), where r c is expressed 
by (98) and the coefficients q\ , q 2 and q% are yielded by 

qi (N e ,x B ) = (l +1 ) [ Ne ^ , (100) 



q 2 (N e ,x L 



(N e .n B ) 7r — arccos N e .n E 
\N K x n B \ \N P x n B \ 



2 l + N K .n r 



-( 1+ ^i-A&? ^ 



q 3 (N e ,X B )^^IlAL 3{l + Ne . nB) + 2(l + 7 f\±^ 
1 - N e .n B { 1 - N e .n J 

tt — arccos N e .n R 



-(l + 7)« 



1 + (1 + 2N e .n B ) 



\N K x n F 



(102) 



Proof of proposition 6 Let x A be a point lying on r s (7V e , x B ). It is clear that 
the part of Y s (N ei x B ) joining x A and x B coincides with a quasi-Minkowskian null 
geodesic path T s (x A ,x B ). So, the impact parameters of T s (N e , x B ) and T s (x A ,x B ) 
are equal. As a consequence, the coefficients qi,q 2 and q% can be obtained as functions 
of N e and x B by taking the limit of equations (94)-(96) when x A recedes towards the 
source of the light ray at infinity, i.e. when r A — >• oo, n A — > — 7V e and N AB — *• N e . 
Taking into account that arccosn^.n B — >• n — arccos N e .n B when n A —J- —N e , we get 
the system of equations (100)-(102). QED. 

The expression found for q$ is new, whereas the expressions obtained for q\ and 
q 2 coincide with previous results found in [5] . 

13. Conclusion 

This paper is devoted to the study of the time transfer function T and the impact 
parameter b corresponding to a photon travelling along a quasi-Minkowskian light ray 
in a static, spherically symmetric spacetime. The main results are the following: 

a) The system of equations (25)- (27) enabling, at least in principle, to determine 
the perturbation terms T^ n > involved in the expansion of the time transfer function 
in power series of G. Such a system appreciably simplifies the approach developed 
in [13] for static, spherically symmetric spacetimes. 

b) The demonstration of the analyticity of the functions T^ when the emission 
and reception points are neither coincident, nor located in diametrically opposite 
directions. 

c) The replacement of the recurrence method outlined in proposition 1 by an 
iterative procedure for solving the integro-differential equation satisfied by the time 
transfer function. This procedure presents the great advantage that only elementary 
integrations which can be performed with any symbolic computer program are required 
whatever the order of approximation. The legitimacy of this approach essentially rests 
on the property of analyticity of the functions T'-™-' . 

d) The explicit calculation of the time transfer function and the impact parameter 
up to the third-order terms in G. The new results brought by equations (87) and (96) 
illustrate the efficiency of the procedure. The expressions obtained for the impact 
parameter up to the third order are extended to the case of a light ray emitted 
at infinity in a given direction. It may be added that the ability to recover the 
expressions of T*- 1 ' and T^ 2 ' found in previous works confirms the reliability of the 
method developed in this paper. 
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It is worth noting that in accordance with equations (40) and (41) in [3], our 
formulae would enable us to derive the triples characterizing the direction of a light 
ray at its points of emission and reception up to the third order. To conclude, it may 
be pointed out that the method could be extended to light rays propagating in the 
equatorial plane of an axisymmctric, rotating body (a Kerr spacetime, e.g.). 



Appendix A. Calculation of T^(x A , x B ) 

To begin with, it may be noted that (78) can be rewritten in the form 
In* = J*. - ( 1 + ^ K + __^£ 4. ( 1 + 7)[(l + 7) 2 + <?iP 

3 3 2 



r: 



r*s 



rs 

\2 , „2 



2s 5 



[3(1 + 7 ) 2 + g 2 ]gi r c r (l+7)tfr; 



2~2 



2s 5 
It is easily seen that 

''' ;: r K3_ _ (1 +j)k nqir c 

r 2 s s 3 rs 3 



(A.1) 



'Vi 



dr 



1 



1 



\Xr 



X, 



l + /i 



K3 - (1 +7) K 



arccos \x 



(A.2) 



The calculation of the three other integrals is more lengthy. For the sake of brevity, 
it is convenient to put 



R 



P = r A r E 



Taking into account (85 a) and (856), and then using relations as follow 
(r B - r A/ u) 3 - (r B n - r A ) 3 = (1 - /j,)(r A + r B )[R 2 (l + f i + M 2 ) 

-p(l + 2 l i- f j?-2 f j?)], 
rl(r B - r A ^f - r 3 A (r B ^i - r A f = R 2 [R 4 + 3pR 2 ^ - 3p 2 (l - /i 2 )], 



3i =(1 + 7) 



2 R 2 + 2p{l + Li) l-/i 



R 2 1 + /J,' 

we are led to 

f*r(l + 7 )[(l + 7 )2 +g 2] r 2 [3(1 + 7) 2 + g 2 ]gircr (l +7 ) g 2 r21 



(A.3) 



(A.4) 
(A.5) 

(A.6) 



r A 



2.S» 



1 



1 



1 



2s 5 

\X R X i 



dr 



(l+ 7 ) 3 / 



ir A r B \r A ' r B J (1 + m) 2 [( r B - r A /u,)(r B n - r A )} 3 ' 
where I is given by the lengthy expression 
/ = R 2 [R 2 +p(l - fi 2 )][R 2 (l +H + M 2 ) - p(l + 2/i - y 2 - 2fi 3 )} 

-(1 - y 2 )[R 2 {2 + n)+p(l - M 2 )][i? 4 + ipR 2 ^ - 3p 2 (l - M 2 )] 

+ (l-/ i ) 2 [i? 2 + 2p(l + M )] 

x [i? 4 (l + 2/i) - pi? 2 (l - 3m - 4/i 2 ) - 3p 2 (l + /u - [i 2 - (i 3 )}. 
Taking into account the relation 

(r B -r A n)(r B (j,-r A ) = R 2 [i-p{l -/i 2 ), 
it is easily checked that 

I 



(A.7) 



(A.8) 
(A.9) 

(A.10) 



\{r B -r A fj,)(r B /j,-r A )] 3 
a result which proves to be spectacularly simple in spite of the apparent complexity 
of the right-hand side of (A.8). Equations (A.2), (A.7) and (A.10) yield (87). 
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